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Siano and Egger Reply: In our recent Letter 
we reported numerically exact path-integral Monte Carlo 
(PIMC) simulation results for the Josephson current 
through a quantum dot in the Kondo regime. Two of 
the central findings of Ref. 0] are as follows. (1) Previ- 
ous studies have missed a minimum in the critical cur- 
rent as a function of A/Tk, where = A u = A is 
the BCS gap in the leads and Tr the Kondo tempera- 
ture for normal leads, and (2) the numerical renormaliza- 
tion group (NRG) method is inaccurate in the crossover 
regime (A/Tr « 1). These findings have been criticized 
by Choi, Lee, Kang and Belzig (CLKB) @. First, they 
point out a mistake in our definition of Tr. As discussed 
below, this error in Ref. |l| does only affect numerical 
values for the various transition points but neither the 
physical picture nor the conclusions drawn in Ref. 0. 
Their second point is that there are strong finite-T ef- 
fects. Here one should first note that PIMC is inher- 
ently a finite-T technique but has the benefit of exact- 
ness. While we did not claim true ground-state results 
in Ref. 0, finite-T predictions are relevant for experi- 
ments, where T w lOOmK corresponding to our value 
for realistic applications. Moreover, our finite-T results 
reproduce the analytical T = solution for A/Tr 3> I 
[lj, in contrast to the assertions of CLKB. Of course, 
at T/A = 0.1, significant thermal smearing in the /(</>) 
curves is present, but this is obvious and also seen in our 
PIMC simulations. Interestingly, a direct comparison of 
the finite-T NRG results shown in Fig. 1 of Ref. to 
exact PIMC simulation data reveals that NRG is not of 
sufficient accuracy for this problem. Problems of NRG 
for the case Al = A^ under discussion here have been 
mentioned also in Ref. 0. 

As pointed out by CLKB, by letting L -> 2r in Eq. (1) 
of Ref. Jj (and only there!), the correct Tr emerges. This 
change implies that equations (8) to (11) in Ref. jlj now 
read: 


plicit comparison of exact PIMC results to the finite-T 


(a/t^v 

(a/TkW 
(a/TkV* 

(A/T^) min 


0.51 ±0.01, 
0.875 ±0.005, 
1.105 ±0.005, 
0.92 ±0.03. 


Furthermore, in Fig. 2 of Ref. 0, we have A/T K « 2.5, 
while in Fig. 3 the corresponding values were A/Tr = 
0.53,0.68, and 0.85. Likewise, in Fig. 4 the correct val- 
ues should read A/T K = 0.85,1.03,1.13, and 11.6. The 
corrected Fig. 5 of Ref. [l| is shown here in Fig. ^ where 
we also added new PIMC data obtained at a lower tem- 
perature in the inset. Fortunately, these changes in the 
numerical values of the various transition points do not 
affect our main conclusions. In particular, Fig. ^ shows 
that the minimum in I c as a function of A/Tr persists 
upon lowering temperature, and is not necessarily wiped 
out by thermal fluctuations. In Fig. [21 we show an ex- 
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FIG. 1: Plot of Fig. 5 of Ref. Q with corrected A/T K units. 
In the inset, new PIMC data at lower T have been added. 
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FIG. 2: Direct comparison between finite-T NRG and PIMC 
results at U/A = 14 and T/A — 2.82. Vertical bars (or the 
symbol size) indicate standard MC error estimates. 


NRG data of CLKB, using precisely the same microscopic 
parameters as in Ref. |2| . Figure|2]indicates that the well- 
established high accuracy of NRG in treating "normal" 
Kondo problems is lost due to BCS gaps in the leads . 
In the crossover regime, A/Tr w 1, NRG is therefore 
unreliable. We conclude that, apart from the mistake re- 
garding Tr posing no fundamental problem, the criticism 
of CLKB is unfounded. 

F. Siano and R. Egger 
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